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00 Abstract 

o 

For type I seesaw and in the basis where the charged lepton and heavy right- 
handed neutrino mass matrices are real and diagonal, four has been shown to be the 
^ ' maximum number of zeros allowed in the neutrino Yukawa coupling matrix yjy. These 

- four zero textures have been classified into two distinct categories. We investigate 

certain phenomenological consequences of these textures within a supersymmetric 
framework. This is done by using conditions implied on elements of the neutrino 
Major ana mass matrix for textures of each category in Y,^. These conditions turn 
out to be stable under radiative corrections. Including the effective mass, which 
appears in neutrinoless double beta decay, along with the usual neutrino masses, 
mixing angles and phases, it is shown analytically and through scatter plots how 
restricted regions in the seesaw parameter space are selected by these conditions. 
We also make consequential statements on the yet unobserved radiative lepton flavor 
violating decays such as /i — > 67. All these decay amplitudes are proportional to 
the moduli of entries of the neutrino Majorana mass matrix. We also show under 
which conditions the low energy CP violation, showing up in neutrino oscillations, is 
directly linked to the CP violation required for producing successful flavor dependent 
and flavor independent lepton asymmetries during leptogenesis. 



1 Introduction 

Neutrino mixing, leptogenesis and (in a supersymmetric framework) radiative lepton fla- 
vor violating decays {ii -^ ^jl), ^ being a charged lepton and i,j being generation indices, 
have generally been acknowledged [1] as important tools to constrain parts of the seesaw [2] 
parameter space. Any study of these tools gets considerably facilitated by the assumption 
of texture zeros being present in the Yukawa coupling matrix Yi, = mc/vu [3-6]. Here Vu 
is the vev of the up-type Higgs and ttt-d is the Dirac neutrino mass matrix. Texture zeros 
also help in relating [3-6] CP violation at low energies to that required for leptogenesis. 
Given the observed pattern of neutrino mixing and assuming no neutrino to be exactly 
massless, four is now known [7] to be the maximum number of zeros allowed in mu within 
the type I seesaw framework. More zero entries in Y^ lead [7] to at least one completely 
unmixed neutrino. This statement is made, of course, in the standard weak basis where the 
right-handed neutrino and charged lepton mass matrices, Mn and m£ respectively, are real 
and diagonal. All such allowed four zero textures in rriD have been completely classified [7] . 
Our aim in this work is to study the implications of these allowed and completely classified 
four zero textures for radiative lepton flavor violating decays as well as for leptogenesis. 
We shall also make observations on related aspects of neutrino mixing and neutrinoless 
double beta decay. 

In type I seesaw [2] the low energy neutrino mass matrix in family space obeys the 'matching 
condition' 

m^ = -ruD M^^ ml = U mf^^ f/^ , (1) 

where U is the PMNS matrix, parameterizable in terms of three angles 612, 6*23, 61^ and 
three phases a,P,6. Thus, 

/ C12C13 S12C13 513 6"*'^ \ 

U = -C23 S12 - S23 si3 C12 e'^ C23 C12 - S23 si3 S12 e'^ S23 ci3 -P , (2) 

V S23 S12 - C23 si3 C12 e*'' -S23 C12 - C23 si3 S12 e*"^ C23C13 / 

where Qj = cos^, Sij = sin 6ij, 6 is the Dirac-type CP- violating phase and the Majorana 
phases a and /3 are contained in the matrix P = diag(l, e*", e**-^^"^-*). Leptonic CP violation, 
that can show up in neutrino oscillation experiments, can be described through a rephasing 
(Jarlskog) invariant quantity given by [6] 

With the parameterization of Eq. ([2]) , one has Jcp = | sin 2^12 sin 2^23 sin26'i3 cos6'i3 sin (5. 
Whereas the presence or the value of any of the phases is currently unknown, the oscillation 
parameters are constrained as follows [8] (see also [9]): 



Ami, = r.67:!:°-(:^°f)xio-^eV^ 



2 , -2.37 ± 0.15 te^) X 10-3 eV^ (inverted ordering) 

^^ ' +2.46 ± 0.15 (Ig:^^) X 10-3 eV^ (normal ordering) , 



sin2 0i2 = 0.32 ±0.02 (+ny , (4) 

„iT,2 /I _ r, /I c- +0.09 /+0.19\ 

sm 6^23 - U.4t)_o.o6 (^-o.iaj ) 

„;„2/i _ r^ n +0.019 /+0.05\ 

sm 6^13 - 0.0 _o 000 [~om) ■ 
The la {3a) ranges around the best-fit values have been given above. 

As in Ref. [7], we consider textures in m^) in the basis in which both mg = diag(me, -m^, rrir) 
and Mr = diag(Mi, M2, M3) are real and diagonal. All fiavor mixing information is thus 
encoded in the Dirac mass matrix mo- The latter can be written in the most general form 
as 



ai e*°i 


02 e*°2 


036*^3 \ 


bi e^^i 


62 e^^2 


63 e*^3 


ci e^^i 


C2 e*^2 


C3e*^3 ) 



mD= \ h e^^i 62 e'^2 ^3 g^/^a . (5) 

\ cie^^i C2e*^2 ^36*^3 ^ 

Here, for each entry, we have listed the real amplitude (a^, bi,Ci) and the corresponding 
phase (aj,/5i,7j) exphcitly. Of course, three of the phases (one per row) can be rephased 
away. This Dirac mass matrix can also be expressed as [10] 



diag 



mD = iU\/mT^R^MR, (6) 

where i? is a complex, orthogonal matrix. This Casas-Ibarra parametrization illustrates 
an important feature: even when the elements of ttIj, and Mr are known, there is still an 
infinite number of Dirac mass matrices leading to the observed low energy phenomenology. 
Other observables need to be used in order to break this degeneracy [11]. 

A well-known strategy to distinguish between different models, leading to the same low 
energy neutrino data, is to make use of Lepton Flavor Violation (LFV) and leptogenesis. 
LEV in supersymmetric seesaw scenarios leads to loop-induced decays such as £j -^ ijj, 
with flavor indices i,j spanning (1 = e, 2 = yU, 3 = r), with the constraint i > j. In 
mSUGRA scenarios, with universal boundary conditions for scalar sparticle mass matrices, 
one obtains the one-loop relation [12] 

BR{ii -^ ij + 7) = (const) BR{ii -^ £j vV) \{moLm^^)ij\^ , (7) 

where the diagonal matrix L is deflned as 

Lui = In 6ki, (S) 

with Mk being the mass of the k^^ right-handed neutrino. The logarithmic coefficient in the 
RHS of Eq. ([8]) takes into account the effect of renormalization group running from a high 



scale Mx to the scale of the respective heavy neutrino masses. The normalization factor 
BR(£j —^ £j vV) in the definition of the branching ratios in Eq. ([7]) is noteworthy. The 
relevant numbers here are BR(r -^ e vV) = 0.178 and BR(r — >• /xz/z/) = 0.174 [13], respec- 
tively. For our later numerical work, we will ignore the small difference between the two. 
We will also take BR(/i -^ e uV) to be unity. Current upper limits on the branching ratios 
for li -^ ij-i are as follows: BR(/i -^ 67) < 1.2 x lO^^^ [14], BR(r -^ 67) < 1.1 x 10"' [15] 
and BR(r -^ fi'-f) < 6.8 x 10^® [16]. One expects these bounds to improve by two to three 
orders of magnitude for BR(/i — > 67) [17] and by one to two orders of magnitude for the 
other branching ratios [18] in the foreseeable future. The unspecified constant in the RHS 
of Eq. ([7]) depends on certain supersymmetry breaking parameters of mSUGRA, specifi- 
cally the universal scalar and gaugino masses and the universal trilinear scalar coupling 
as well as on tan/5. However, we are not interested here in the exact magnitude of the 
branching ratios. We shall instead study the vanishing of certain branching ratios which 
for ii — ^ ^j + 7 turn out to be proportional to the square of the i, jth element of the low 
energy mass matrix m,^. 

In principle, the above analysis could be extended also to other lepton flavor violating 
processes, such as /i-e conversion in nuclei [19]. Current experimental limits on this pro- 
cess, however, are expected to be improved considerably only much after stronger limits 
on £i — > ijj have been made available. If the photon penguin contribution dominates 
the LEV diagrams, as happens for the case under study, a good estimate for the ratio 
of BR(/x — i> 67) to the rate of /i-e conversion is 0{l/a), where a is the electromagnetic 
fine structure constant. In particular, the rate of /x-e conversion is also proportional to 
{niD Lnvj^)i2. Hence, if in one of the scenarios to be discussed BR(/i -^ 67) vanishes, /x-e 
conversion will be absent as well. Note, moreover, that since only one conversion channel 
{n -^ e) is experimentally accesible for the conversion process, no potentially testable dou- 
ble ratios can be given. Eor these reasons, our focus here is on the ii —>■ ijj decays. 

The other important aspect of seesaw phenomenology is leptogenesis. Of particular interest 
are the decay asymmetries [20,21] that depend explicitly on the charged lepton flavor: 



e — 






(9) 



where 



^ = Im 



J 



^D)iai^D)aii^D^D)ii , Ji'' = lm{m^ij)ia{mD)aj{m^D^D)ji- (10) 



.U (^ \ l^\ 



It is evident that Xg = -Z'^^ and J^^ = -Jf^. In the MSSM, the function /(x) has the 
form [22] 

2 , /l + a;N 



f{x) = ^/x \- In ( ) 

L — nr. V T / 



1 — X 

4 



11 



We have given quite general expressions above for the decay asymmetries that can accom- 
modate any nontrivial role played by flavor effects [24] . Thus ef describes the decay of a 
heavy right-handed neutrino of mass Mj into a charged lepton of flavor a = e, /i, r. When 
the lowest-mass heavy neutrino is much lighter than the other two, i.e. Mi <C M2,3, the 
lepton asymmetry is dominated by the decay of this lightest of the heavy neutrinos. In this 
case f{Mj/Mf) ~ —3 Mi/Mj. Moreover, only the first term proportional to X"- in Eq. iQ 
is relevant then since the second term proportional to J7^" is suppressed by an additional 
power of Mi/Mj. Note furthermore that the second term in Eq. (Q vanishes when one 
sums over flavors to obtain the flavor independent decay asymmetry: 



E^ 



.a a 



^[r(iv,^0y-r(iv, ^0t/j- 



1 



1 



/3 



1 



5]lm 



[m^D^ofij 



fiMf/Mf) 



:i2) 



where we have defined 



^71 vl {m^^mD)i 



E^- 



■X,- 



u ' 



(13) 



We note here already that for all 72 four zero textures that we study the J7]"- vanish. The 
expressions given above for the decay asymmetries are valid for the MSSM. Their flavor 
structure is, however, identical to that of just the Standard Model. 



Equally important in leptogenesis are effective mass parameters that are responsible for the 
wash-out. We assume that a single heavy neutrino of mass Mi is relevant for leptogenesis. 
Then every decay asymmetry e" gets washed out by an effective mass 



[m 



Ln 



a {mo) 



al 



Ml 



(14) 



Moreover, the wash-out can be estimated by inserting this effective mass in the function [25] 



ri{x) 



8.25 X 10-3 eV 



+ 



X 



X 



2 X 10-4 eV 



1.16^ 



(15) 



The summation of m'^ over the flavor index a yields mi, which is the relevant parameter 
for the wash-out of Si. The final baryon asymmetry is [23-25] 



Y, 



B 



—0.01 £177(7721) 
-0.003 ({e\ + e^i) r, (||(m^ + <)) + ^I V (fi^I 
^ -0.003 (e\ T] (i§ ml) + £^ r^ (f|| m^) + ^I r^ (f|| ml 



one-flavor , 

two-flavor, (16) 

three-flavor . 



Here we have given separate expressions for one-, two- and three- flavored leptogenesis. 
The three-flavor case occurs for Mi (1 + tan^ (3) < 10^ GeV, the one-flavor case for Mi (1 + 
tan^ (3) > 10^^ GeV, and the two-flavor case (with the tau-flavor decoupling first and the 
sum of electron- and muon-flavors, which act indistinguishably) applies in between. 

LFV and leptogenesis provide means of breaking degeneracies in the seesaw parameter 
space. This comes about since their dependence on the seesaw parameters is complemen- 
tary to that of m^. A related issue is the question of circumstances under which there is a 
connection between high and low energy CP violation, i.e., between the phases responsible 
for leptogenesis and the ones responsible for CP asymmetries in neutrino oscillations. Inas- 
much as texture zeros simplify this process, the motivation behind the present study is to 
phenomenologically investigate the extent of this degeneracy breaking for all allowed Dirac 
mass matrices with four zero textures classified in [7] . The rest of the paper is organized 
as follows: in Section [2] the two categories, (i) and (ii), of the four zero textures in tjid 
are recapitulated and the radiative stability of the corresponding conditions on m^ is em- 
phasized. In Section [3] we discuss the phenomenology of these textures, focusing on lepton 
mixing and the ratio of ratios in radiative LFV decays as well as on leptogenesis, wash-out 
factors and the basis invariant Jarlskog CP- violating parameter Jcp! subsections 13.11 and 
13.21 cover categories (i) and (ii) respectively. The final Section H] contains a summary of our 
results and conclusions derived therefrom. 

2 The Two Categories of Four Zero Textures 

It will be helpful to provide first a summary of the classification of the four zero textures. 
As enumerated in Ref . [7] , there are 72 allowed textures of this kind: 

(i) 54 textures in which two rows of rriD are orthogonal element by element. They can 
be further divided into three subclasses containing 18 matrices each: 

(ia) 18 textures in which the first and second row are orthogonal element by element, 
resulting in 

{m^)i2 = {rn^)2i = 0, (17) 

i.e., the vanishing of the off-diagonal 12 (or efi) entry of the effective neutrino 
Majorana mass matrix; 

(ib) 18 textures in which the first and third row are orthogonal element by element, 
resulting in the vanishing off-diagonal element condition 

("^^)i3 = ("^1^)31 = 0, (18) 

i.e., the vanishing of the off-diagonal 13 (or er) entry of the effective neutrino 
Majorana mass matrix; 



(ic) 18 textures in which the second and third row are orthogonal element by element, 
resulting in 

{m^)23 = {m^)32 = , (19) 

i.e., the vanishing of the off-diagonal 23 (or fir) entry of the effective neutrino 
Majorana mass matrix; 

(ii) 18 textures in which two columns of itld are orthogonal element by element. They 
can be further divided into three subclasses containing 6 matrices each: 

(iia) 6 textures with two zeros in the first row, resulting in the vanishing sub- 
determinant conditions 

|(m^)ii (m^)23| - \{m^)2i {rn^)i3\ = arg{(m,,)ii (m,,)23 {rn^T2i {m^Yn} = ; 

(20) 

(iib) 6 textures with two zeros in the second row, resulting in the vanishing sub- 
determinant conditions 

|(m^)22 (m^)i3| - 1(1^1^)12 ("^i.)23| = arg{(m^)22 {m^)i3 (m^Tu ('^1^)23} = ! 

(21) 

(iic) 6 textures with two zeros in the third row, resulting in the vanishing sub- 
determinant conditions 

|(m^)33 (m^)i2| - |(m^)i3 {171^)32] = arg {(772^)33 (m^)i2 (m^)*3 (m^)*^ = . 

(22) 

We now wish to comment on the question of the dependence of such results, as presented 
above, on the energy scale. In general, elements of the matrix Y^, change with energy in a 
coupled way due to radiative corrections leading to renormalization group (RG) running. 
Thus an element, which vanishes at low energies, can certainly develop a significant nonzero 
value at a very high energy. Our postulate is that four zeros are present in Y^, at energies 
relevant to oscillation experiments. Within a reasonable accuracy, such can also be taken 
to be the case then for radiative LFV decays. This assumption is, however, generally not 
valid for leptogenesis which we take to operate at Mi > 10^ GeV. In a grand unified theory, 
of course, Y^ would originate at Mx — 10^^ GeV and would need to be evolved down to an 
energy scale below the Z-mass rriz- If texture zeros are due to some (yet unknown) flavor 
symmetry at some high scale, one would need to assume that the said zeros are protected 
by the same symmetry during the RG running. If such is the case, our statements would 
continue to hold without modification. 

Let us nevertheless point out a particularly interesting feature of the consequences for 
the neutrino Majorana mass matrix of the textures in 772^5 under consideration. These 
conditions on m,y are stable under RG running. If one performs the running from the high 



scale Mx to the low scale mz-, then for mi, this can be taken into account by multiplying 
each matrix element {m^)ij by a factor [27]. The latter is given by (1 + ej) (1 + e^), where 

m? . M 



1 "^ JC 

m ■ 



with rTi3,2,i = fnT,fi,e being the charged lepton masses. The parameter c is given by 3/2 in 
the SM and by —(1 + tan^/5) in the MSSM. The multiplicative nature of this correction 
ensures that a vanishing element of rriu stays vanishing. Thus the consequence for rui, 
of every texture under category (i) is safe under RG running. It is, additionally, rather 
surprising that the corresponding consequences for rrii, from all textures in category (ii) 
are also unharmed by the RG running from radiative corrections. For instance, consider 
category (iia) and perform the corrections in condition (120|) . The multiplication of {my)ij 
with the factors (1 + e^) (1 + e^) leads to 

|(m^)ii (mj,)23| - \{mu)2i {mu)i?\ -^ |("^i^)ii {my)'2z\ - l("^i^)2i {^y)'i?,\ 

= |(m^)ii (m^)23| (1 + ei)2 (1 + ea) (1 + £3) - |(m^)2i (m^)i3| (1 + ei) (1 + £2) (1 + ei) (1 + £3) 

= (1 + eif (1 + £2) (1 + eg) (l(m^)ii (m^)23| - \{m^)2i (m^)i3|) , 

so that condition fl2Ul) remains unchanged. The result is identical for categories (iib) and 
(iic), the matrix indices 1, 2, 3 appearing the same number of times on both sides of the 
conditions (|2T1) and (l22l) . This defines an interesting class of "RG invariants" . 

3 The Phenomenology of Four Zero Textures 

3.1 Category (i) 

We first discuss some general issues concerning the phenomenology of the textures under 
consideration. To start with, take the subclass in category (i) in which rows i and j (7^ i) 
of rriD are orthogonal, element by element: 

{m^)ij = {m^)ji = . (23) 

It also follows that for this subclass of four zero textures, 

{mo K m^j,)ij = , (24) 

where K is any diagonal matrix. The immediate implication is that the branching ratio 
for the decay ii — > ij'j is zero for these textures: 

BR(£i -^ Ij-f) = . (25) 

Therefore all such textures of mo would be excluded by any future experimental observa- 
tion of this decay modqj. 



^The requirement of a vanishing (to_d to|))i2 can lead via 2-loop effects to a lower limit on BR(yLt — > 67), 
connected to the product of the branching ratios of r — > /i7 and r — > 67 [26] . Obviously the 2-loop induced 
branching ratio is very small. 
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Another interesting feature is that, in the textures of category (i) and for an arbitrary 
j ^ i-, the relation 

\{mDL7ll}J:))ij\(x\{m^,)ij\ (26) 

holds. However, it is not possible to construct predictive ratios of branching ratios from 
this, unless the heavy neutrino masses are known. Consider as an example the texture 



rriD 



/ 02 as e 
6i 

V ci C2 e^T2 



las 



which belongs to category (ia) and leads to {rny)i2 = BR(yU — > 67) 
branching ratios for the decays r — ^ 67 and t —>■ fij are governed by 



(27) 



0. The non-zero 



{moLrn^jj] 



13 



On Co Ln and 



''2 "^2 



{mDLm^^)23 =hlclL 



'2 

'1 5 



(28) 



respectively. The low energy Majorana mass matrix is 



/ 



m^ 



2 „2ia3 



M2 



+ 



a, e 



M. 





Ml 



a2C2e 



272 



M2 
61 Ci 



M- 



Mi 



+ 



cse 



2172 



(29) 



M2 J 



It follows that 



BR(r^e7) |(m,)i3|' /M2L2^' 



BR(r^/i7) |(m,)23|' KM.LJ ' ^^°^ 

Without any further information about the heavy neutrino masses, one is unable to predict 
this ratio. The same feature is valid for all textures in category (i). We will not work out 
here all 72 possibilities. 



Leptogenesis - either of the unfiavored or of the flavor dependent variety - is quite possible 
in general for such textures. As mentioned earlier, we work under the assumption that 
Ml <^ M2,3 so that one needs to consider only the decay of A'^i. Thus, among the parameters 
responsible for leptogenesis, i is always 1. Generalization to the more general situation, 
including M2,3 is, however, straightforward. It turns out that for the textures under study 
here all J'{j vanish and we are left only with the T^. Specific textures in this subclass will 
of course have some of the X^ zero in case the appropriate elements of rriD happen to lead 
to this. In the example from Eq. fl27|) the only non-zero decay asymmetry is 



T 



12 



cl c\ sin 272 . 



(31) 



It is interesting to ask under what circumstances the "leptogenesis phase" 72 is responsible 
for low energy leptonic CP violation as well. Evaluating the invariant in Eq. ([3]), which 
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Category (ia) 

inverted ordering 



Category (ia) 

normal orderinc 




•> 4'-* 












0.0001 



0.001 



0.01 



m, [eV] 
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Figure 1: Category (ia) or {m^)i2 = 0: scatter plots of the effective mass versus tlie smallest 
mass in case of an inverted mass ordering (left) and of sin^ 26'i3 versus the smallest mass 
for a normal mass ordering (right). For the first plot we have also given the general upper 
and lower limit of the effective mass when the currently allowed 3a values of the oscillation 
parameters are used. The corresponding plots for category (ib) look basically identical. 



describes CP violation in neutrino oscillations, gives 

— CL2 Oi Ci C2 



Am2i ^mL AmL Jqp 



(cl Ml sin 2a3 + (bj + cf) M2 sin 2(03 - 72 



Mf Ml M3 

- {{bl + cl)al + blcl) Ms sm2^2 



(32) 

It follows that the leptogenesis phase is related to the low energy Dirac phase when the 
conditions 



bi + ci]a^2 + K4jM3 >a: 



cl Ml sin 20:3 



(bl + cl)M2 sin 2{a3 



72, 



(33) 



are fulfilled. Similar considerations can be made for all other textures in category (i). Note 
that in the basis that we are working, the Dirac mass matrix from the example in Eq. ( 127|) 
contains only two physical phases, which is one less than the number of low energy phases 
in m^. This facilitates the connection of low and high energy CP violation. We stress here 
that all four zero textures contain two physical phases. For this reason all 72 candidates 
have - in analogy to the conditions in Eq. (155]) - the possibility [6] that the low energy 
leptonic CP violating phases are the ones responsible for leptogenesis. 
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Figure 2: Category (ic) or (m;y)23 
case of a normal mass ordering. 



Category (ic) 

nonnal ordering 










ii [eV] 



0: scatter plot of sin^ 6'i3 versus the smallest mass in 



Let us next discuss the neutrino mixing properties of category (i), in which case the low 
energy neutrino mass matrix contains a vanishing off-diagonal entry. The phenomenology 
of mass matrices rrii, having one single texture zero was analyzed in Ref. [28]. Consider 
first the case of (mj,)i2 = 0, i.e., category (ia). The first thing to note is that zero f/es is 
incompatible with the condition of vanishing (?7ij,)i2 [28]. To elaborate, expanding in terms 
of It/esl one finds 



{'mu)i2 = {rn2e 



2ia 



rrii 



COS ^12 sin 6*12 cos ^23 
'^^"-COS^^ 



(34) 



+e'^ sin ^23 {e''" m^ - sin" Ou m2 e'*" - cos" 6^ mi j \Ue3\ 
plus higher order terms of If/esI- The magnitude of the bracketed part of the zeroth order 



term is bounded from below roughly by JAtjIq for a normal hierarchy, by ^AttIq/ JAm\ 
for an inverted hierarchy and by ^AmQ/rriQ for quasi-degenerate neutrinos with an average 
mass rriQ. The 12-element therefore cannot vanish for [/e3 = 0, but can vanish if Ues departs 
from zero. It is clear from the above expression that in case of an inverted hierarchy and 
for quasi-degenerate neutrinos the bracketed part of the zeroth order term in Eq. flMl) has 
to be small in order to allow the term of order 6'i3 to cancel it. This in turn means that 
sin a has to be close to zero. In this case there are almost no cancellations in the effective 
mass (m) = |(m,^)ii| governing neutrinoless double beta decay and we have 



(35) 



(m) ~ J Am\ cos^ ^13 inverted hierarchy , 
mo cos^ 6*13 < (m) < mo quasi- degeneracy . 

The above value and range have to be compared with the lower limits, arising from max- 



imal effects of Major ana phases, i.e.. 



(m)^ 



cos^i3 cos2^i2 A/Am^ E'^iid (m) 
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008^6*13 008 26*12^0, respectively. The left panel of Fig. [T] shows for category (ia) a scatter 
plot of the effective mass versus the smallest mass in case of an inverted mass ordering. It 
is clearly seen that the largest possible values of (m) are mostly populated. We also show- 
in the right panel a scatter plot of sin^ 26*13 versus the smallest mass for a normal mass 
ordering. In order to generate this and other plots to be presented later, we have varied 
all neutrino parameters within their allowed Sa-ranges quoted in Eq. (jl]). The results for 
category (ib), i.e., (mj,)i3 = 0, are basically identical to the ones for category (ia) [28]. 
Formally one can obtain the 13-entry of rrii, from the 12-entry by replacing in the latter 
sin 6*23 with cos 6*23 and cos ^23 with —sin 6*23. 



Category (ic) 

normal orderina 



Category (ic) 

inverted ordering 




nij [eV] 



mj [eV] 



Figure 3: Category (ic) or (7/1,^)23 = 0: scatter plots of the effective mass versus the smallest 
in case of a normal mass ordering (left) and an inverted ordering (right). We have also 
given the general upper and lower limits of the effective mass when the currently allowed 
3cr values of the oscillation parameters are used. 



Turning to category (ic), it can be shown that, for 6*13 = and a normal mass ordering, the 
23-element of m,y cannot vanish unless neutrino masses are above several eV [28]. Fig. [2] 
illustrates this by showing the distribution of the smallest mass against sin^ 2^13 in case of 
a normal mass ordering. For quasi-degenerate neutrinos, we can express the mass matrix 
element as 



l"^^J23 



mo cos 6*23 sin 6*23 I (e *° cos 6*12 + sin 9 



12 



,2i(/3+<5) 



where for simplicity we have also set ^13 to zero. It is clear that, in order to make 
("^i/)23 vanish, the expression in the brackets should be very small, or that the relations 
sin a ~ sin(/5 + (5) ^ should hold. This leads again to little cancellation in the effective 
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mass, as is obvious from Fig. [31 

We close this Section by commenting on the possibihty that more than one entry of the 
low energy mass matrix vanishes. In this respect it is known that, in the mass diagonal 
charged lepton basis, two is the maximum number of vanishing elements allowed in the 
neutrino mass matrix [29]. We refer to Ref. [29] for details on the phenomenology of these 
cases. It is sufficient to note here that seven of those cases exist, namely the simultaneous 
vanishing of the 11- and 12-, the 11- and 13-, the 12- and 22-, the 13- and 22-, the 13- 
and 33-, the 12- and 33-, and finally the 22- and 33-entries. In category (ia), in which 
('^i/)i2 = 0, there is the possibility that in addition (m,^)ii, (7/1,^)22 or {171^)33 can be zero. 
This in turn means that the branching ratios of the decays r — ;> 67 and r — > /i7, which 
depend on |(mi,)i3p and |(m^)23p, respectively, are guaranteed to be non-zero. The same 
is true for category (ib), in which case {m^)i3 = 0, and for which again only {m^)ii, {my)22 
or {my)33 can be zero. BR(/i -^ 67) and BR(r -^ fi'j) then are non-zero because there are 
proportional to the non-zero |(m,^)i2p and \{m^)23\^, respectively. In case of category (ic), 
or {mu)23 = 0, it turns out that no other mass matrix element can vanish, and therefore 
BR(/i —>■ 67) and BR(r — > 67) are necessarily non-zero. 

3.2 Category (ii) 

Now, let us similarly consider the subclass of category (ii) in which columns / and k (7^ /) 
are orthogonal, element by element. For this subclass of four zero textures, the relation 

(m^mD)/fc = (36) 

applies. This means that the terms XjJ, and ei vanish in the respective fiavored and unfia- 
vored heavy neutrino decay asymmetries, cf. Eqs. (jHl) and (IT^ . 

It will be illuminating to explicitly see the different seesaw induced physical effects that 
appear in a given subclass of textures. Tabled] has been made for this purpose by consid- 
ering category (iia). The leftmost column contains the Dirac mass matrices for particular 
textures of this category after one phase per row has been rotated away. This is followed in 
subsequent columns by expressions for only those heavy neutrino decay asymmetries which 
are nonzero and contribute accordingly to leptogenesis as well as their corresponding wash- 
out factors. For simplicity, we work as before under the assumption that Mi -C M2^3 so 
that one needs to consider only the decay of Ni, i.e., i is always 1. However, this is an 
inessential assumption. Our tables can easily be generalized to the cases with j 7^ i. Inter- 
estingly, we find again that the JT^" are all zero. The entries in the last column of Table [1] 
are expressions describing the invariant Jcp relevant to CP violation that can be observed 
in neutrino oscillations. The latter is always the sum of three terms and we have chosen 
not to list common proportionality factors. Again, in analogy to the discussion in the 
previous Subsection, a comparison of the entries in the leptogenesis and the Jcp columns 
is instructive. It shows how there could be a one-to-one correspondence between the lep- 
togenesis phase and the low energy Dirac phase S. For instance, consider the first row in 
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Table [T], for which only el contributes to leptogenesis, and for which the CP phase denoted 
by 7i is crucial. The same phase can control low energy (Dirac) CP violation provided the 
condition 



bj Ml M3 > |ci M2 M3 sin 2/32 1 



+ bl + cl)MiM2 sin2(/52 



7iJ 



(37) 



is fulfilled. In the last two rows of Tabled] there are two nonzero decay asymmetries, which 
lead to more possibilities. In this respect, one may note that hierarchical heavy neutrinos 
lead to a suppression of X^g with respect to 1^2 by a factor M2/M3. Tables [2] and [3] show the 
leptogenesis related phenomenology of categories (iib) and (iic), respectively. The matrices 
can be obtained from the ones of category (iia) by interchanging the first with the second 
(first with the third) row. 

Even though none of the branching ratios is guaranteed to vanish in this class of textures, 
there exists an interesting feature. For all 18 matrices in category (ii) the following relation 
applies for i 7^ j: 

\{mDLm^jj)ij\ oc \im^)ij\ , (38) 

i.e., the branching ratios for the decays ii — >■ ij'y are proportional to the square of the 
modulus of the ij element of the low energy mass matrix. This is similar to the situation 
in category (i). The crucial difference is however that in category (ii) the ratios of branching 
ratios are related to ratios of neutrino mass matrix elements without any ambiguity coming 
from the unknown values of the heavy Majorana neutrino masses. Consider again the 
example from the first row of Table [H The low energy Majorana mass matrix is given as 



/ 



m,. 



03 

Ms 



M3 



Mo 



+ 



ft2 

M, 



a3C3 

M3 
&3C3 

Ms 

c2e2ni 



\ 



v 



Ml 



+ 



(39) 



Ms ; 



The relevant expressions for LFV are found to be 



'm£i Lnir 



D)12 



%^3-^3 ' 



[mDLm^ij)i3 



ttg Cs Lg and 



{mnLm^: 



23 



63 Cs Ls 



Therefore, the ratios of branching ratios are unambiguously given by the ratios of the 
corresponding mass matrix elements. This "minimal lepton flavor violation" scenario in 
principle allows one to predict the rates from measurable low energy mass matrix elements. 
Indeed, from Eq. ([7]), we see that the following relations now hold in category (iia): 



BR(r 



BR(r 
BR(/x -^ 67) 
BR(r -> 67) 



^^) ^, 


("^i.)23 


2 


^7) ~ 
1 


5 

('Ti,.)l2 


BR(r- 


-> euu) 


im^)i3 



(40) 
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Figure 4: Category (iia): scatter plots of sin^ 26'i3 
sin^ 26^12 versus sin^ 6*13 (for mi = 0.001 eV) in case 
the conditions |(m^)ii {171^)231 = I("^t/)i2 ("T'i/)i3| and 
are fulfilled. 



versus the smallest mass mi and of 
of a normal mass ordering and when 
Im{(m,,)ii (mj,)23 (mj,)i2 {rn^Yis} = 



Eq. (140|) is valid for all textures of category (ii) . 

The three sub-categories in category (ii) have in addition correlations between the low 
energy mass matrix elements, given in Eqs. (12UI f2T\ [2^ . which lead to 



and 
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Figure 5: Category (iia): scatter plot of the effective mass (tti) versus the smallest 
mass ?7i3 in case of an inverted mass ordering and when the conditions \{vti^x\ [j^v)2'i\ = 
li'm,,) 12 {rriu) isl and Im{(m,y)ii (my)23 (m,y)*2 (^!/)i3} = are fulfilled. The solid (green) 
lines correspond to the upper and lower limit of the effective mass when the currently 
allowed 3a ranges of the oscillation parameters are used. 



Because the neutrino mass matrix obeys /x-r (or 2-3) symmetry to a good approximation, 
\{'mu)i2/{m^)i3\'^ is typically 1. Therefore, BR(r -^ 67) ~ 0.178 BR(/x -^ 67). With the 
current upper limit of 1.2 x 10~^^ on BR(/i — > 67), and an expected improvement of at 
most two orders of magnitude on the limit of BR(r -^ 67) < 1.1 x 10^^, it follows that in 
this scenario r — ^ 67 will probably not be observed in a foreseeable future. Since it turns 
out that BR(r — > fi'j) ~ BR(r — > 67), the same is true for the decay r -^ /i7. To discuss 
further the phenomenology of category (iia), we show in Fig. H] scatter plots for a normal 
mass ordering in case the correlation fl20l) holds. Obtaining analytical correlations for this 
example is very complicated. If we set mi = ^13 = and 623 = 7r/4, then Eq. fl20|) leads to 
I jAm\/AmQ = cos^ 6*12 cos2(a — P — 5) < 1, which is incompatible with the available 
data. This implies the presence of a lower limit on sin^2^i3, which indeed turns out to be 



'13, 



roughly 0.04 (0.06) for mi = 0.001 (0) eV. This lower limit vanishes when the neutrino 
mass increases, which can clearly be seen in Fig. HI 



With regard to an inverted mass hierarchy, we find no correlation for the neutrino oscillation 
parameters. Instead, the effective mass (m) = |(wii,)ii| governing neutrinoless double beta 
decay is constrained by the correlation in Eq. fl20l) . Note that in the definition of the 
correlation for category (iia) the effective mass ((wi^)ii) appears explicitly. The influence 
on (m) can be estimated by considering the equality |(my)ii (777.^)23^ = \{rn^)2i ("^i^)i3p 



in an approximate manner. Neglecting 777.3, setting 777.2 



777i 



A777^ 



A; 



and inserting 
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sin^^23 = |) "W'e obtain from it the condition 
(Amir 



1 - 2 sin^ 2^12 sin^ a) + 0{\Ue3n = . 



I 



Hence, sin^ a ~ 1/(8 sin^ 9i2 cos^ 612) (which is equal to ^ if sin^ 6*12 = |) and inserting 
this in the effective mass leads to 

cos^ 61 



W^^^^V^^- (43) 

This has to be compared with the general lower and upper limits on (m), which are 



cos 26*12 cos^^i3 JAm\ and cos^ 6*13 JAm\, respectively. Fig. [5] shows how the effective 
mass as a function of the smallest mass m^ has considerably less spread than without the 
correlation Eq. (!20|) . For the other two conditions in Eqs. (12T!) and (122|) we did not find 
any interesting correlations between the neutrino observables. 

Turning to the ratios of branching ratios, Eqs. fj4Tl) and fj42l) are phenomenologically very 
interesting. In Fig. [6] we show as a function of sin^6'23, the ratio of the branching ratios 
of r ^ 67 and r ^ /i7 in the left panel and oi fi ^ ej and r ^ 67 in the right panel. 
These numbers are equal to |(7Ti,^)i3p/|(m,y)23p and (|(mj,)i2p/|(m,^)i3p)/BR(r — >• eVu), 
respectively. Note that, for normal ordering, the /x-r block of rrii, is usually larger than the 
elements of the e-row, while for the inverted ordering case, all elements of rrii, are of similar 
magnitude. This explains why in the left panel the predictions for the inverted ordering are 
higher compared to that for the normal ordering. In addition, since BR(r -^ eVv) = 0.178, 
the ratios shown in the right panel are larger than those in the left panel. To be more 
specific, let us consider a simplified example. For inverted ordering, assuming 6*13 = m^ = 0, 
1712 — 'mi and sin^ 6^12 = I, we obtain 



\{mu)i?. 
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|(m,,)23p COs2 6'23 

which explains the mild increase of this ratio as a function of sin^ 6'23 in the left panel of 
Fig. [6l On the other hand, the ratio BR(/i -^ e7)/BR(r -^ 67) is seen to be proportional 
to cot^ ^23 for the same set of assumptions. This therefore results in the decrease seen for 
inverted ordering in the right panel of Fig. [61 For normal ordering and with the same set of 
assumptions one can show that |(?7^i/)i2p/|("^i/)i3p oc tan^6'23. This explains the increase 
of BR(/i -^ e7)/BR(r — > 67) with sin^ ^23 in the right panel of Fig. O Note that, for the 
inverted mass ordering BR(r — > e7)~BR(r — * 7x7), while BR(/i -^ e7)~BR(r — > 67) is 
true irrespective of the neutrino mass spectrum. Therefore, we do not expect BR(r — > 67) 
to be observed in the forthcoming experiments [18]. Similarly for the inverted ordering, the 
predicted BR(r — > /i7) is not expected to be checked experimentally in the next generation 
experiments [18]. The scatter plots of BR(/i — > e7)/BR(r — > 67) in all sub-categories of 
category (ii) look very similar and the other cases show no interesting correlations. 
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Figure 6: Category (iia): scatter plots of BR(r — * e7)/BR(r -^ fi'y) (left) and 
BR(yU — > 67)/ BR(r — > 67) (right) against sin^2^23 when the conditions |(mi.)ii (r?7,,^)23| = 
\{^TT'u)i2 ("^i/)i3| and Im{(m,^)ii {m^)23 ("^1^)12 ("^1^)13} = are fulfilled. The red circles are 
for the normal ordering, the green squares for the inverted one. 



Finally, we discuss the possibility and implications of vanishing branching ratios, and hence 
the low energy mass matrix elements. Since the basic criteria for category (ii) comprise 
of zero sub-determinant conditions (cf. Eq. (I20|) - (l22|) ). one can easily note that in this 
case, the vanishing of one element of m^ would necessitate the vanishing of another mass 
matrix element as well. Therefore, for category (ii) we cannot have just one zero m^ 
texture, and hence just one vanishing branching ratio. The allowed two zero textures 
of m„ have been extensively studied in the literature [29]. One can check that none of 
the phenomenologically viable two zero textures allow (mi,)23 to be zero. An immediate 
consequence of this is that for category (ii), BR(r — > fi'y) ^ 0. Among the sub-categories, 
we note that for category (iia), if BR(/i — > 67) = (implying {m^)i2 = 0), then {m^)n = 0. 
Therefore, for this case a vanishing BR(/i — > 67) predicts vanishing neutrinoless double beta 
decay - this in turn is possibile only for normal mass ordering with mi ~ 0.005 eV. For 
category (iib), the only allowed two zero texture is with {171,^)12 = and {my)22 = 0. For 
this category a vanishing {my)i^ would make {rny)i2 = and this case is strongly disfavored 
by the data. Therefore, BR(r -^ 67) 7^ for this case while BR(/i -^ ej) could go to zero. 
On the other hand, for category (iic), the {m^)i2 7^ condition is imposed by the data, 
implying that BR(/i — > 67) 7^ 0, while BR(r -^ 67) could go to zero if both (mjy)i3 = 
and {171^)22 = simultaneously. 
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4 Summary and Conclusions 

We have looked at phenomenological constraints from the two maximally allowed cate- 
gories of four zero textures in Fj, in the basis in which the mass matrices M/j and rrii are 
diagonal. Our framework is that of the supersymmetric type I seesaw and we have exam- 
ined the consequences on seesaw parameters of the conditions imposed on elements of the 
neutrino Majorana mass matrix rrii, by either category of textures. We have included the 
effective Majorana mass (m), appearing in neutrinoless nuclear double beta decay, in the 
list of seesaw parameters studied. Our use of those conditions is reliable in that the latter 
are radiatively stable, being invariant under RG running. 

We have demonstrated via various scatter plots how restricted regions in the seesaw pa- 
rameter space are selected by the said conditions. For the textures of each category, we 
have further derived a number of results on radiative LFV decays, several of them with 
observable consequences. For instance, any observation of the decay ii -^ ij'-f would rule 
out those category (i) textures which imply that {m^)ij = {m^)ji vanishes. In general, 
all 72 four zero textures predict that the branching ratios of ii -^ ij'y are proportional to 
the absolute values squared of the mass matrix element {'m^)ij times a function of heavy 
Majorana neutrino masses. For category (ii) this function is the same for all i,j and the 
ratios of branching ratios are directly given by ratios of low energy mass matrix elements. 
For category (ii), we have been able to generate sample scatter plots for ratios of branching 
ratios such as BR(/i — ^ e7)/BR(r — > ej) and BR(r -^ e7)/BR(r -^ fi'y) against sin^^23- 
We have also obtained results with physical consequences for leptogenesis from these tex- 
tures. An example is that, for both categories (i) and (ii), we have derived conditions 
which fix the functional form of the leptogenesis phase purely in terms of the low energy 
Dirac phase without invoking any other phase from among the seesaw parameters. This 
is always possible for the four zero textures because their number of physical phases is 
two. We have also provided additional information by tabulating quantities relevant to 
lepton flavor asymmetries, wash-out factors and the Jarlskog invariant for each texture of 
category (ii). 

In conclusion, we have highlighted the rich phenomenological structure of the allowed four 
zero textures in Yi, deflned in the charged lepton and right-handed neutrino mass diagonal 
basis. The consequent conditions on rriiy are radiatively stable. These have been shown to 
lead to signiflcant reductions in the seesaw parameter space. 
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